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I.  Introduction 


In  November  1982  MIT  proposed  to  conduct  a  definition  study  of  a 
coordinated  group  of  experiments  to  measure  the  electromagnetic  and  plasma 
environment  of  a  space  shuttle  in  polar  orbit.  This  group  of  experiments  was 

intended  to  provide  scientific  measurements  and  data  sufficient  to  interpret 

% 

the  results  of  engineering  experiments  carried  on  the  Interaction  Measurement 
Payload  for  Shuttle  (IMPS).  The  original  objectives  of  the  study  were  as 
fol  1  ows : 

1.  Coordinate,  with  AFGL,  the  definition  of  the  scientific  and 
programmatic  requirements  of  the  investigation,  including  an 
understanding  of  any  active  plasma  and  particle  experiments 
envisioned  in  the  future. 

2.  Review  current  literature  covering  theory,  observations  and 
instrumentation,  including  visits  to  appropriate  universities  and 
other  organizations. 

3.  Obtain  a  thorough  understanding  of  the  Space  Trail  portation  System, 
its  operation,  structure,  matenials  and  contaminants,  with  an 
emphasis  on  the  effects  the  Shuttle  will  have  on  the  environment  and 
vice  versa. 

4.  Develop  preliminary  psuedo-stati onary  models  of  the  Shuttle's  plasma 
and  electromagnetic  environment  as  it  would  be  influenced  by  the 
environment  of  a  typical  polar  orbit. 

5.  Explore  the  perturbations  of  the  environment  as  it  is  affected  by  the 
dynamical  interactions  with  external  sources  (solar  and/or 
geomagnetic  events)  or  internal  sources  (active  plasma  experiments, 
debris  dumps,  terminator  crossings). 

6.  Recomnend  a  package  of  coordinated  experiments  to  be  flown  on  an 
early  polar  orbiting  Shuttle  to  investigate  the  Shuttle's  local 
envi ronment . 

7.  Present  an  implementation  plan  to  procure  these  experiments, 
integrate  them  into  the  Shuttle,  operate  them  in  orbit,  collect  data, 
and  finally  analyze,  interpret  and  publish  the  results  in  a  timely 
fashion. 

8.  Provide  a  budgetary  estimate  of  item  No.  7.  above  based  upon  an 
expected  flight  of  the  Envi ronmental  Monitoring  Package  in  FY1987. 
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The  original  proposal  was  funded  for  the  period  February  18,  1983  through 


November  30,  1983  and  later  extended  through  September  30,  1984.  Program 
delays,  reorientation  and  rescoping  of  the  IMPS  program  and  an  inability  on 
the  part  of  MIT  to  find  qualified  personnel  to  work  on  the  program  resulted  in 
some  changes  in  emphasis  and  a  reduction  in  scope  of  the  original  objectives. 
During  the  early  phase  of  the  project,  theoretical  work  was  started  on  a 
better  understanding  and  general  theoretical  treatment  of  the  dispersion 
relations  for  a  cold  plasma  and  on  a  general  theoretical  treatment  of  the 
interaction  of  a  moving  conductor  with  an  ambient  magnetized  plasma 
(objectives  3,  4,  and  5).  A  major  effort  was  devoted  to  objective  6  (the 
definition  of  a  coordinated  set  of  scientific  support  experiments. 

1 1 .  Recommended  Scientific  Monitoring  Payload 

A  report  and  description  of  the  recommended  support  experiments  was 
completed  and  included  in  a  progress  report  submitted  to  AFGL  in  April  1984. 
Obviously,  the  appropriate  payload  of  scientific  support  experiments  depends 
on  the  engineering  experiments  finally  selected  for  IMPS  and  some  modifica¬ 
tions  to  the  April  report  are  included  in  the  final  report  of  work  to  evaluate 
tht'  engineering  experiments  proposed  by  JPL  which  was  carried  out  by  MIT  under 
AFGL  contract  F19628-83-K -0033. 

A  major  point  stressed  in  the  MIT  payload  recommendations  was  that  the 
scientific  support  experiments  required  measurements  not  just  in  the  payload 
bay  but  as  a  function  of  distance  from  the  STS.  We  suggested  initially  that 
the  Spartan  free  flyer  developed  at  GSFC  be  used  for  this  purpose  and  this 
recommendation  led  directly  to  the  idea  of  using  the  German  SPASS  vehicle  to 
carry  out  the  scientific  support  measurements.  Note  that  in  the  JPL  "Report 
on  Final  Recommendations  for  IMPS  Engineering/ Science  Payload"  and  in  the  MIT 


Final  Report  on  the  evaluation  of  the  JPL  recommendations,  the  ensemble  of 
scientific  support  experiments  is  referred  to  as  EIM,  the  Environmental  and 
Interactions  Monitor. 

III.  Theoretical  Studies 

One  of  the  most  important  questions  concerns  the  interaction  of  the  STS 
and  its  payloads  with  the  ambient  magnetized  plasma.  In  general  a  conducting 
body  (even  an  insulated  one)  moving  relative  to  a  magnetized  plasma  generates 
a  variety  of  plasma  waves  which  not  only  are  effects  of  scientific  interest 
but  may  also  result  in  practical  engineering  consequences  and  limitations.  A 
related  problem  is  the  propagation  of  a  given  wave  mode  in  the  plasma  as  a 
function  of  physical  parameters. 

Obviously,  radiation  of  electromagnetic  signals  by  the  STS  or  other 
conducting  bodies  in  orbit  is  of  direct  interest  for  the  IMPS  program  for  many 
reasons,  e.g.,  noise  generation,  satellite  drag,  charging  effects,  etc.  For 
this  reason  Barnett  and  Olbert  undertook  a  general  study  of  the  problem  which 
they  describe  in  Appendix  11.  This  work  has  been  support''d  first  by  AFGL 
under  contract  F19628-83-R-0030  and  in  its  final  stages  b'  NASA/JPL  under  JPL 
contract  963733.  A  final  paper  describing  this  research  las  been  completed 
and  submitted  to  the  Journal  of  Geophysical  Research. 

A  complete  analytic  and  numerical  treatment  of  the  phase  and  group 
velocities  of  waves  in  the  cold  plasma  approximation  was  undertaken  by  Dr. 
Zdzislaw  Musielak.  This  work  was  completed  in  October  1984  and  is  available 
as  a  Center  for  Space  Research  Technical  Report,  CSR-TR-84-3.  This  TR 
constitutes  Appendix  III  ot  the  present  final  report. 


IV. 


Future  Studies 


To  fully  understand  the  interaction  of  the  STS  and  other  large  structures 
witn  the  ambient  medium,  the  theoretical  work  described  above  must  be  expanded 
to  include  major  modifications  of  the  local  STS  plasma  environment  which 
result  from  large  quantities  of  neutral  gas  emitted  by  the  STS  and  cold  plasma 
dispersion  relations  need  to  be  modified  to  include  warm  and  hot  plasmas 
actually  present  during  typical  STS  polar  orbits.  Eventually  it  should  be 
possible  to  predict  the  near  and  far  radiation  fields  expected  and  to  design 
realistic  experiments  to  measure  the  actual  plasma  waves  generated  by  motion 
of  the  STS  and  other  large  structures  through  the  ambient  plasma. 
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THE  ENVIRONMENT  MONITORING  PAYLOAD 


I.  INTRODUCTION 

The  purpose  of  the  Environment  Monitoring  Payload  (EMP)  is  to  determine 
the  in  situ  plasaia  and  electromagnetic  environment  in  and  around  the  Shuttle 
in  polar  orbit  and  to  measure  the  interaction  of  the  Shuttle  with  that 
environment. 

The  rear-Earth  plasma  and  electromagnetic  environment  has  been  studied 
and  monitored  for  well  over  two  decades.  The  physics  of  the  ionosphere  and  of 
the  ionospheric-magnetospheric  interaction  are  at  present  fairly  well 
understood.  For  example,  known  mechanisms  for  production  and  loss  of  charged 
and  neutral  particles  lead  to  models  of  the  low  and  raid-  atitude  ionosphere 
which  are  in  reasonable  agreement  with  observations.  At  higher  latitudes, 
extensiive  spacecraft  measurements  have  defined  both  tne  energy  spectrum  and 
the  raagnetospheric  source  regions  for  those  precipitati nt  particles  which 
produce  the  visible  aurora,  and  thus  at  least  the  qualititive  physics  of  the 
auroral  oval  is  in  hand. 

With  the  advent  of  the  Space  Transportation  System  (i.e.,  the  Shuttle)  a 
different  aspect  of  the  near-Earth's  environment  must  be  analyzed  and 
monitored  since,  for  the  first  time,  the  size  of  the  space  vehicle  structure 
and  its  own  atmosphere  will  substantially  influence  and  possibly  dominate  the 
in  situ  plasma  and  electromagnetic  environment.  Some  preliminary  measurements 
have  been  and  will  be  made  on  early  STS  flights  from  the  Kennedy  Space  Center 
(KSC).  However,  a  Shuttle  in  polar  orbit  will  experience  a  far  more  varied 
environment  due  to  the  much  wider  range  of  geomagnetic  latitudes  (e.g.,  the 
STS  will  sample  the  high  fluxes  of  energetic  precipitating  particles 


associated  with  the  auroral  zone).  The  interaction  of  the  STS  with  this 
environment  is  largely  unknown. 

II.  BRIEF  DESCRIPTION  OF  THE  IN  SITU  ENVIRONMENT  AND  POSSIBLE  SHUTTLE 
INTERACTIONS _ 

The  lowest  practical  altitude  for  a  polar  orbiting  Shuttle  is  about  200 

kilometers.  At  this  altitude  and  at  raid- latitudes,  the  Shuttle  will  be  in  the 

FI  region  of  the  Ionosphere.  The  various  geophysical  regions  of  interest  that 

such  a  Shuttle  will  sample  are;  1.)  the  polar  cap;  2.)  the  auroral  zone;  and 

3. )  the  low  to  middle  geomagnetic  latitudes  below  the  auroral  zone. 

In  the  polar  cap  region,  O'*^  is  probably  the  dominant  ion,  and  the  ambient 

tnermal  plasma  has  a  density  of  10^  to  10^  electrons/ cc,  with  temperatures  of 

a  few  tenths  of  an  eV.  There  are  soft  precipitating  electrons  in  this  region 

(with  energies  from  10  eV  to  100  eV)  which  are  thought  to  be  magnetosheath 

electrons  which  have  gained  entry  into  the  mantle  region  and  are  subsequently 

precipitated.  In  addition,  there  is  a  plasma  outflow  (the  polar  wind)  of  0^. 

In  tie  nightside  auroral  zone,  the  ambient  plasma  densities  are  an  order  of 

magnitude  higher  than  in  the  polar  cap  region.  This  increased  density  results 

from  ionization  of  neutrals  by  the  high  fluxes  of  hard  precipitating  particles 

in  the  1  to  10  keV  range  (predominately  electrons  but  with  proton  fluxes  at 

lOi  of  the  electron  levels).  In  the  worst  cases,  auroral  fluxes  can  reach 
10  2 

2  X  10  electrons/ cm  sec  steradian.  There  is  also  an  appreciable  flux  of 

electrons  backscattered  up  the  magnetic  field  lines,  and  copious  numbers  of 

secondary  electrons  are  produced  in  the  10  eV  to  100  eV  range.  On  the 

dayside,  photoionization  produces  electrons  with  10  eV  to  50  eV  energies,  with 

higher  ambient  densities  compared  to  the  nightside,  and  there  is  a  softer 

spectrum  of  precipitating  particles.  Finally,  in  the  low  to  mid-latitudes, 

there  is  no  high  energy  electron  component,  and  the  electron  densities  (all 

5  6 

due  to  photoionization)  run  from  10  to  10  electrons/cc  in  the  dayside. 


The  interaction  of  the  Shuttle  with  the  plasma  environment  in  a  low  Earth 
orbit  is  complicated.  Theoretically,  there  is  some  understanding  of 
steady-state  sa  tellite/plasma  interactions  at,  for  example,  geostationary 
orbit.  Here,  the  Debye  length  of  the  plasma,  and  the  electron  and  proton 
gyroradii  of  the  tb.ermal  particles,  are  large  compared  to  the  dimensions  of  a 
small  satellite.  In  addition,  orbital  velocities  are  much  less  than  typical 
ion  thermal  speeds.  Such  ambient  plasma  conditions  at  geostationary  orbit 
allow  one  to  make  a  number  of  approximations  in  modeling  spacecraft  charging, 
for  example,  that  the  plasma  is  unmagnetized,  and  that  the  satellite  is  at 
rest  with  respect  to  the  plasma.  In  low-Earth  orbit,  in  constrast, 
temperatures  are  much  lower  and  magnetic  field  strengths  and  densities  are 
much  larger.  As  a  consequence,  the  Debye  length  is  only  millimeters  to 
centimeters,  the  gyroradii  of  the  ambient  electrons  are  centimeters,  and 
gyroradii  of  typical  ions  are  meters.  In  addition,  tiic  Hhut*le  orbital  speed 
i 'ti  km/' sec)  greatly  exceeds  the  ambient  or  ion  thermal  ;;peeds  (~1  km/sec). 
Thus,  in  modeling  spacecraft  charging,  one  cannot  assume  an  unmagnetized 
plasma,  and  one  cannot  ignore  wake  effects.  Because  of  these  complications, 
there  are  major  conceptual  problems  with  any  existing  quasi-static  charging 
code  as  applied  to  the  Shuttle  in  low-Earth  orbit.  In  addition  the  worst  case 
in  terms  of  the  most  severe  charging  environment  for  the  Shuttle  will 
’undoubtedly  be  auroral  currents,  which  have  spatial  scales  down  to  tens  of 
meters.  At  Shuttle  orbital  speeds,  the  Shuttle  may  saranle  such  features  r.o 
rapidly  that  quasi-static  approximations  to  charging  nay  not  be  valid.  In 
this  and  other  areas,  there  are  major  theoretical  problems  in  modeling  the  STS 
plasma  interaction,  and  such  theoretical  developments  will  undoubtedly  have  to 
be  guided  and  verified  by  experimental  results. 


In  addition  to  the  question  of  the  interaction  of  the  Shuttle  with  the 
ambient  plasma,  there  is  the  question  of  changes  in  ambient  conditions 
introduced  by  the  presence  of  the  Shuttle,  e.g. ,  contamination  due  to  venting, 
outgassing  or  thruster  firings,  or  >;harging  of  the  Shuttle  due  to  "active" 
experiments  such  as  electron  beam  experiments. 

Recent  flights  of  the  Space  Shuttle  have  shown  that  the  local  environment 
is  dominated  by  contamination  from  ;he  Arbiter  itself.  This  contamination  is 
in  the  form  of  gases,  particulates,  and  glows  that  are  either  generated  by  the 
Shuttle-atmosphere  interaction  or  a  result  of  integration  procedures  and 
pre-launch  conditions.  These  topics  have  been  discussed  in  detail  at  a 
variety  of  meetings  (AIAA  meeting,  51  October  1983  to  2  November  1983, 
Washington,  D.C.;  AGU  meeting,  5  December  1983,  San  Francisco;  URSI  meeting, 

11  January  1984  to  13  January  1984,  Boulder),  in  journal  articles  (e.g.,  see 
Geophys.  Res.  Lett.,  V10,  ffll,  1983  for  Shuttle  glow  studies),  and  in 
technical  reports  (e.g.,  contamination  reports  in  NASA-82-457  and 
NASA-82-524).  The  following  limited  summary  highlights  areas  of  concern  for 
the  Shuttle  environment. 

-4  -7 

.  local  gas  ggessures  vary  from  10  to  10  Torr  compared  to  the 
ambient  10  Torr, 

.  water  vapor  (H^O)  is  the  principal  contaminant  but  high  levels  of 
helium  (He)  are  also  noted, 

.  high  levels  of  electrostatic  plasma  waves,  energetic  ions,  and 

suprathermal  electrons  have  been  measured — perhaps  related  to  Shuttle 
glow, 

particulate  levels  are  significant  primarily  during  periods  of  ascent 
and  descent  and  then  only  in  the  <  5  micron  range. 

Naturally,  many  of  the  above  measurements  are  greatly  influenced  by  and 
dependent  on  Shuttle  attitude  (e.g.,  bay-to-sun  versus  passive  thermal 
con'rol),  orbital  maneuvers  (e.g.,  thruster  firings),  and  operations  (e.g.. 


10 


I  sf 


’  n  PI",  • 


water  dumps).  The  number  of  varying  parameters  complicates  the  analysis  and 
the  effects  of  the  Shuttle- induced  atmosphere  on  sensitive  experimental 
measurements  (e.g.,  infared  radiation  or  ambient  plasma  measurements)  are  far 
from  being  understood. 

All  NASA  Shuttle  flights  have  been  limitod  to  the  low  and  mid  latitudes 
and  have  not  traversed  the  auroral  zones  and  polar  caps.  The  auroral  region 
is  characterized  by  intense  auroral  electron  precipitation.  DMSP  satellites 
within  this  region  have  shown  that  such  fluxes  of  charged  particles  can  lead 
to  large  electric  potentials  on  the  surface  ol  the  vehicle.  Further, 
laboratory  measurements  of  typical  space-qualified  materials,  which  have  been 
irradiated  with  auroral-type  electors,  have  shown  that  significant  charge 
deposition  can  occur  within  exposed  dielectrics.  Both  situations  are 
hazardous  since  the  possibility  of  arcing  exists,  and  this  arcing  is  known  to 
upset  on-board  electronics.  In  the  case  of  the  Shuttle,  charge  effects  are 
expected  to  be  especially  acute  since  the  variety  and  nature  of  the  surface 
materials  are  susceptible  to  differential  and  deep-dielectric  charging. 

III.  RELATIONSHIP  TO  OTHER  PROGRAMS 

The  thrust  of  the  EMP  is  to  measure  the  interaction  of  the  Shuttle  and 
its  contaminants  with  the  local  plasma  and  electromagnetic  environment  in  a 
polar  orbit.  Existing  data  sets  and  scheduled  NASA  projects  do  not 
specifically  address  the  auroral  environment.  Specifically,  the  OSS-1 
Pathfiinder  Mission  probed  the  near-Shuttle  environment  but  at  low  latitudes. 
SPACELAB-II  will  have  an  increased  measuring  capability  but  is  still  limited 
in  latitude.  The  instrumentation  proposed  for  EMP  is  spec  .flcally  geared 
towards  measurements  in  the  auroral  zones. 
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The  earliest  possible  flight  opportunities  are  reconnnended  t'or  EMP  in 
order  to  use  the  information  gathered  in  the  design  and  flight  planning  of 
sensitive  Air  Force  Systems.  Without  knowledge  of  the  contamination  levels 
and  effects  of  auroral  particle  precipitation.  Air  Force  systems  may  be 
neeclessly  impacted.  It  is,  therefore,  strongly  recommended  that  these 
experiments  be  conducted  ASAP  from  a  polar-orbiting  Shuttle  mission. 

With  the  exception  of  EMP,  no  DoD  payloads  have  been  specifically  defined 
to  study  the  contamination  levels  and  the  Shuttle-auroral  environment 
interaction  on  a  polar-orbiting  Shuttle  mission.  However,  some  other  programs 
are  complimentary  and  their  results  should  be  incorporated  into  the  overall 
investigative  program. 


IMPS  (AFGL) 


SC I ED  (AFGL) 


SEPS  (NRL) 


PIE-II  (ONR) 


is  geared  towards  the  engineering  aspects  or 
consequences  of  the  environment.  It  is  planned, 
however,  that  IMPS  will  use  the  EMP  to  accomplish  its 
mission  objectives  by  measuring  the  probable  causes  of 
those  consequences. 

is  concerned  with  large  structures  within  the  space 
environment.  The  auroral  environment  is  not  the 
principal  concern  here  and  the  instrument  compliment 
reflects  this.  The  experimental  techniques  to  be  used 
by  SCIED  to  measure  the  wake/ram  effects  of  large  bodies 
in  space  are  different  from  EMP  but  will  be  useful  in 
understanding  the  interaction. 

is  concerned  with  plasma  irregularities  and  their  effect 
on  communication  and  surveillance  systems.  No  provision 
is  made  for  detailed  auroral  investigations. 

has  provisions  to  study  the  effects  of  the  auroral 
environment  on  spacecraft  but  is  an  active  experiment 
attempting  to  determine  the  effects  of  man-induced 
plasma  processes  rather  than  those  of  natural  processes. 


IV.  THE  EMP  PROGRAM 

To  fully  explore  the  interaction  of  the  Shuttle  with  the  polar- 
magnetospheric  environment  will  require  a  free-flying  probe  capable  of 
maneuvering  throughout  the  three-dimensional  space  in  the  vicinity  of  the 
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Shuttle.  Other  dimensions  (time,  orbital  parameters,  shuttle  aspect  angles, 
e.g.,  ram/wake,  etc.)  must  be  varied  in  order  that  a  full  characterization  of 
the  interaction  be  obtained. 

It  is  recognized  that  many  practical  reasons  will  restrict  the  program 
from  achieving  all  of  these  spatial  and  temporal  measurements  with  a  full 
compliment  of  experiments.  It  is  equally  important  to  realize  that  even  with 
a  full  complement  of  measurements  it  is  still  a  very  difficult  non-linear 
theoretical  problem  to  fit  a  model  to  the  results.  The  program  must  begin 
somewhere,  and  as  long  as  we  keep  our  sights  on  the  ultimate  payload  and  do 
not  expect  ultimate  results  from  a  leaser  payload,  then  a  somewhat  less 
capable  payload  can  be  initiated  to  obtain  some  early  results  in  the  1987  time 
frame. 

A.  The  Early,  Captive  Program;  SACS 

A  group  of  Shuttle  Auroral  and  Contamination  Sensors  (SACS)  can  be 
assembled  on  a  pallet-like  structure  to  measure  the  gaseous,  particulate, 
and  plasma  densities  and  the  electromagnetic  fields  surrounding  the  Shuttle. 
While  these  sensors  will  not  give  a  true  and  accurate  measure  of  the 
interaction,  even  if  deployed  by  the  remote  manipulating  system  (RMS),  they 
will  nevertheless  give  a  fair  indication  of  the  in  situ  conditions  in  the 
Shuttle's  cargo  bay,  particularly  in  regard  to  local  contaminants.  The 
Shuttle  should  present  as  many  different  aspect  angles  with  respect  to  its 
orbital  velocity  as  the  time-line  can  accommodate  so  that  ram/wake  effects  can 
be  studied  at  various  locations  in  the  polar  orbit. 

The  following  complement  of  instruments  should  be  included  in  the  SACS 
package: 

Mass  Spectrometer  -  The  ion/neutral  quadrupole  mass  spectrometer  obtains 
composition  measurements  of  all  incoming  gases.  Then  by  taking  a  second  set 
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of  measurements  with  a  retarding  voltage  on  a  grid,  one  is  able  to 
discriminate  against  those  gases  that  are  being  carried  along  with  the 
Shuttle.  The  difference  between  the  two  measurements  give  an  estimate  of  the 
Shuttle  contamination.  The  preferred  direction  is  into  the  RAM.  The 
f le Id-of-v iew  is  -  20°  and  tne  mass  range  is  1  to  64  AMU. 

Energetic  Particle  Detectors  -  This  package  includes  an  electrostatic 
analyzer  and  two  solid-state  detectors  which  will  measure  the  energy  spectra 
of  electrons  and  ions  in  the  ranges  of  1  eV  to  150  keV  and  1  eV  to  1  MeV, 
respectively.  Such  instrumentation  is  commonly  used  to  detect  auroral 
particle  fluxes  and  to  detenine  spacecraft  charging.  These  particle 
detectors  are  optimized  when  viewing  in  the  zenith  direction  although  the 
nadir  direction  allows  measurements  of  atmospherical ly-backscattered  particles. 

UV/Visible  Spectrometer  -  This  spectrometer  consists  of  two  telescopes 
(FL  =  230  mm)  and  two  1 /8-meter  Ebert-Fastie  scanning  spectrometers  mounted  on 
a  scan  platform  to  change  thi  look  angle.  A  latching  npchanism  protects  the 
gears  during  launch  and  landing  and  a  retractable  cover  at  the  telescope 
entrance  is  used  to  protect  the  optics  and  provide  an  internal  reflecting 
surface  for  the  liveness  lamp.  To  prevent  detector  damage  by  intense  solar 
radiation,  there  is  a  sun  sensor  which,  when  activated,  turns  off  the  high 
voltage  and  closes  the  shutter.  The  wavelength  range  is  determined  by  an 
optical  grating  and  a  photon  counting  photomultiplier  tube. 

Stereo  Cameras  -  Two  l6-mni  Photo-Sonics  IPI  cameras  are  to  be  mounted  in 
a  stereo  configuration  along  with  an  optional  strobe  light.  Each  camera  has  a 
self-contained  400-ft  film  magazine  containing  high-speed,  high-contrast  film. 
The  cameras  have  internal  heaters  and  are  enclosed  in  a  pressurized  container. 
IRIG-B  timing,  generated  by  the  Orbiter,  is  recorded  on  each  film  frame. 
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Pressure  Gauge  -  This  measurement  technique  is  a  gas-dischar<?e  manometer 


using  a  self-sustained,  low-pressure  discharge  created  by  means  of  crossed 
magnetic  and  electric  fields.  The  electric  current  extracted  from  the 
discharge  region  is  proportional  to  the  sensor  internal  pressure  which,  in 
turn,  is  established  by  ambient  gas  flow  from  the  environment  into  the  sensor 
volume.  Laboratory  calibrations  and  kinetic  gas  theory  provide  the  basis  for 
the  interpretation  of  sensor  output  in  terms  of  ambient  gas  flow  or  density. 

Near-IR  Photometer  -  A  PbS  detector  assembly  with  mechanical  chopper  and 

circular  variable  filter  between  2  and  3  microns  observes  near-IR  emissions 

from  contaminants.  Laboratory  calibrations  establish  the  relation  among 

2 

voltage  output,  filter  position,  and  the  observed  photon  flux  (watts/ cm  -ster). 
The  photometer  has  a  t  20°  f ield-of-view. 

Plasma  Probes  -  1.)  The  Langmuir  probe  measures  the  local  plasma  density 
and  temperature.  As  presently  defined,  it  sweeps  from  -4V  to  ♦AV  in  4  seconds 
and  has  j^edestal  offsets  of  0,  4,  6,  12,  and  16V.  In  'ne  fixed-potential  mode 
the  probe  measures  variations  in  electron  density. 

2.)  Either  2  or  5  orthogonal  sets  of  dipole  antennas 
will  measure  the  local  DC  and  AC  electric  fields  within  the  Shuttle  bay.  Two 
swept-frequency  receivers  operating  in  the  ranges  0  to  60  KHz  and  0  to  5  MHz 
cover  all  expected  plasma  wave  emissions  excited  by  the  Shuttle-atmospheric 
interactions. 

Magnetic  Fields  -  1.)  The  fluxgate  magnetometer  measures  the  full  DC  and 
near-DC  (<40  Hz)  magnetic  fields  in  the  vicinity  of  the  Orbiter.  The 
sensitivity  of  this  instrument  is  sufficient  to  detect  auroral  field-aligned 
currentr  but  this  measurement  may  be  heavily  impacted  by  Shuttle-generated 


noise . 


2. )  A  3-axls  search  coil  is  used  to  measure  the  AC 
magnetic  field  in  the  same  frequency  range  as  the  AC  electric  system.  It  may 
be  possible  to  use  the  same  swept-f requency  receiver  as  the  electric  field  by 
appropriate  switching  of  the  input  to  the  receiver. 

Cameras  -  Within  the  aft  flight  deck  conventional  TV  and  still  cameras 
are  to  be  used  to  photograph  the  aurora  during  geomagnetically  active  periods. 
A  wide-angle  camera  within  the  bay  is  to  be  used  to  photograph  the  aurora 
under  a  more  controlled  situation.  Clearly,  for  this  latter  case  the  Orbiter 
bay  should  be  nadir  pointing.  A  simple  grating  may  be  used  in  front  of  the 
camera  to  yield  spectral  information. 
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Table  1 


Experiment  Complement  and  Requirements 


Dimensions  Dimensions  Weight 


Item 

Stowed  (cm)  Deployed  (cm) 

(kg) 

Ejected 

Recovery 

Mass  Spectrometer 

N/A 

15x7x7  bxbx6 

12.7 

N/A 

N/A 

Energetic  Particle  Det 

N/A 

25x25x30 

14.2 

N/A 

N/A 

UV/VIS  Spectrometer 

N/A 

11x16x21 

18.6 

N/A 

N/A 

Stereo  Cameras 

N/A 

TBD 

33.6 

N/A 

N/A 

Pressure  Cauge 

N/A 

6x5x5  7x4x5 

4.1 

N/A 

N/A 

Near-IR  Photometer 

N/A 

4x6x4 

2.3 

N/A 

N/A 

Plasma  Probes 

N/A 

7x6x7  14x4x4 

5.0 

'  N/A 

N/A 

Magnetic  Fields 

N/A 

7x6x7 

5.0 

N/A 

N/A 

Cameras 

N/A 

10x10x10 

4.5 

N/A 

N/A 

Nominal  Power 

150  - 

200  watts 

Peak 

Power 

180  - 

250  watts 

Energy 

6  - 

10  kWh 

Total 

.  Mass 

100  - 

150  kg 

Real 

Time  Data  Rate 

13  - 

30  kbps 

Data 

Storage 

10'° 

bits 

Payload  Specialist  Duties; 

.  Turning  power  ON/OFF  and  monitoring  critical  functions. 

.  During  magnetically  active  peiiods  when  there  are  intense  auroral 
features,  the  astronauts  may  l«  requested  to  photograph  the  displays 
using  on-board  cameras  (still  and  TV),  To  this  end,  it  may  be 
necessary  to  liave  the  astronauts  change  film  canisters  in  addition  to 
doing  the  photography. 


B.  The  Free  Floater  EMP 

The  only  meaningful  way  of  determining  the  interaction  of  the  Shuttle 
with  the  plasma  and  electromagnetic  environment  is  to  make  measurements 
throughout  the  three-dimensional  region  in  the  vicinity  of  the  Shuttle.  One 
needs  to  get  far  enough  from  the  Shuttle  to  avoid  its  own  atmosphere  of 
contaminants  and  explore  the  effects  that  the  Shuttle  fas  on  waves,  particles 
and  fields. 

The  complement  of  experiments  is  basically  the  same  as  in  the  SACS 
package.  The  major  difference  is  that  it  now  has  to  be  totally 
self-sufficient  and  recoverable.  There  are  several  possible  candidate 
carriers  but  perhaps  the  closest  and  most  readily  available  is  the  SPARTAN, 
currently  in  production  at  GSFC.  It  has  all  the  desirable  features  of  a 
free-floating,  recoverable.  Shuttle-compatible  carrier,  i.e.,  it  has  an 
approved  safety  package,  structural  analysis,  support  systemr  and  STS 
interface,  including  grapple  and  release  mechanism.  It  is  a  three-axis 
stabilized  carrier  but  lacks  translational  ability.  This  is  not  a  major 
drawback  since  the  Shuttle  could  maneuver  around  it,  if  the  SPARTAN  could  not 
be  so  enhanced. 

The  first  SPARTAN  will  be  released  from  the  Shuttle  this  August  1984  and 
four  others  are  currently  in  production  or  in  the  late  planning  stages.  While 
the  SPARTAN  may  not  be  the  only  available  carrier  (SPAS,  PDP,  Multi-probe  are 
other  possibilities),  it  is  well  conceived  and  a  complete  system  which  will  be 
flight  proven  this  year.  It  should  be  considered  very  seriously  for  the 
free-floating  EMP.  The  full  capabilities  of  SPARTAN  are  described  in  the 
attached  document  from  GSFC. 
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Appendix  II 


Radiation  of  Plasma  Waves  by  a  Conducting  fedy 
Moving  Through  a  Magnetized  Plasma 
by  A.  Barnett  and  S.  Olbert 

We  are  studying  the  interaction  between  a  moving  conducting  body  and  a 
magnetized  plasma.  It  is  well  known  that  a  large  conducting  body  moving 
slowly  ac ross  magnetic  field  lines  radiates  low  frequency  (Alfven)  waves 
(Orell,  Foley  A  Ruderman  (1965)).  In  this  work  we  estimate  the  power 
radiated  by  a  moving  sphere  or  cylinder  of  arbitrary  size  at  all 
frequenc ies  . 

Previous  workers  in  the  field  of  spacecraft  interactions  with  plasmas 
concerned  themselves  primarily  with  computing  the  electrical  charging  of 
spacecraft.  They  used  models  which  neglected  the  magnetic  field  and  the 
"motional”  electric  field.  Instead,  they  concentrated  on  detail ed  model  i ng 
of  the  electrical  properties  and  geometry  of  the  spacecraft.  They  computed 
the  current  balance  by  computing  individual  particle  trajectories  and 
applying  the  Liouville  theorem.  The  magnetic  field  was  neglected  because 
its  inclusion  makes  the  particle  trajectories  prohibitively  difficult  to 
compute.  For  a  small  satellite  in  geosynchronous  orbit,  this  type  of  model 
is  adequate.  For  large  structures  in  low  Earth  orbit,  these  models  neglect 
the  dominant  effect. 

The  first  work  to  concentrate  on  effects  due  to  the  magnetic  field  and 
the  "motional"  electric  field  was  that  of  Drell,  Foley  &  Ruderman.  Using 
"cold  plasma"  theory  in  the  MHO  limit,  they  showed  that  a  large  conducting 
body  in  low  Earth  orbit  radiates  Alfven  waves,  and  proposed  that  the  energy 
loss  due  to  this  mechanism  was  the  cause  of  the  anomalously  fast  rate  of 
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decay  of  the  orbit  of  the  Echo  I  weather  satellite.  Their  explanation  did 
not  gain  wide  acceptance,  however,  and  it  was  not  until  direct  in  situ 
observation  of  the  magnetic  field  perturbation  and  later  of  the  plasma  flow 
perturbation  due  to  the  A1  fven  wave  generated  by  the  satellite  lo  in  the 
Jovian  magnetosphere  that  the  phenomenon  was  first  observed. 

The  present  work  is  basically  an  extension  of  the  DFR  paper.  We 
consider  not  only  the  MHO  regime  but  also  waves  of  higher  frequencies. 

We  formulate  the  problem  in  a  self-consistent  fashion  for  the  passive 
state,  and  we  estimate  the  power  radiated  and  the  total  electric  current  of 
the  system.  We  make  certain  simplifying  assumptions,  similar  to  those  made 
by  Drell  et  al :  1)  The  conducting  body  moves  "slowly"  relative  to  the 
sationary  plasma  ("slow"  compared  to  c^,  the  Al  fven  speed  of  the  ambient 
medium).  2)  The  scale  lengths  of  the  gradients  in  the  unperturbed  plasma 
and  magnetic  field  are  large  compared  to  the  size  of  the  moving  conductor 
so  that,  in  effect,  the  unperturbed  medium  may  be  treated  as  uniform  and 
homogeneous.  3)  The  velocity  of  the  conductor,  V,  remains  constant  and  may 
be  treated  as  a  given  perameter.  Moreover,  ^  lies  in  the  plane 
perpendicular  to  the  ambient  magnetic  field. 

Although  our  work  is  closely  related  to  the  issue  of  "s:)acecraft 
charging"  widely  covered  in  the  literature,  our  approach  is  vastly 
different  from  that  taken  by  the  workers  in  that  field.  We  completely 
bypass  the  discussion  of  the  microscopic  details  near  the  surface  of  the 
conductor  and  concentrate  on  the  macroscpic  collective  effects  relatively 
far  from  the  body. 

Our  formulation  is  based  tpon  the  familiar  procedures  of  plasma 
theory  with  one  important  modification:  we  augnent  the  general  equation  for 
the  electric  field  by  a  "source  term"  involving  the  electric  current  within 
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the  moving  conductor.  More  specifically,  we  start  with  the  following 
equation  in  the  -Fourier  space; 

u  4i:i 

l(x(l(xE)  +  _  K•^  =  -  _  (1) 

2 

c  .J 

where  i(  is  the  wave  vector  ?  is  the  electric  field,  K  is  the  dielectric 
tensor,  is  the  angular  frequency,  c  is  the  ^eed  of  light  in  vacuum,  and 
>5^  is  the  source  current.  Applying  Ohm's  law  in  the  conductor,  we  can 
write  an  integral  equation  which  determines  3^  in  a  self-consistent  manner. 
Since  the  integral  equation  is  difficult  to  solve,  we  proceed  to  study 
solutions  to  Equation  1  for  a  known  source  current,  and  we  infer  the 
spatial  dependence  of  the  source  current  fron  physical  argirnents. 

We  begin  by  studying  the  dielectric  tensor  and  the  dispersion 
relation,  including  constraints  on  them  imposed  by  our  problem.  Ke  seek 
solutions  which  (a)  do  not  depend  on  time  v^en  viewed  from  a  frame  of 
reference  attached  to  the  moving  conductor  and  ( b)  involve  only  real 
propagation  vectors  k  in  the  Fourier  space.  The  latter  constraint  is 
directly  coupled  to  the  primary  purpose  of  the  work,  namely,  calculation  of 
the  energy  loss  of  a  moving  conducting  body  solely  due  to  the  radiation  of 
plasma  waves.  We  find  that  both  of  the  above  conditions  can  be  satisfied 
only  by  that  mode  of  plasma  waves  that  has  nearly  vanishing  phase  velocity 
for  some  specific  direction  and  magnitude  of  the  propagation  vector  It.  (Ely 
"nearly  vanishing"  we  mean  a  phase  velocity  smaller  than  the  velocity  of 
the  conductor.)  The  entire  treatment  is  very  reminiscent  of  the  classical 
treatment  of  Cherenkov  radiation  produced  by  fast  moving  bodies  in  an 
ordinary  diectric  medium.  Because  of  the  enormous  complexity  of  the 
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diectric  tensor  of  plasmas,  we  restrict  our  explicit  evaluations  to  the 
cold  plasma  approximation.  We  believe  that  this  approximation  is 
sufficient  for  exploratory  purposes.  We  do  not  aim  at  precision  but  rather 
to  point  out  the  existence  of  new  phenomena . 

We  find,  in  the  cold  plasma  approximation,  that  waves  which  satisfy 
the  above  stated  conditions  exist  only  in  the  following  three  frequency 
bands ; 

0  <  u  <  2a 

w,  L,  <  w  <  2b 

Ln  6 

<  u)<  2C 

where  ,  Q^,  and  Up  are,  respectively,  the  ion  and  electron 

cyclotron  frequencies,  the  lower  and  upper  hybrid  frequencies  and  the 
plasma  frequency.  The  first  band  consists  of  the  A1  fven  waves,  already 
studied  by  Orel!  et  al .  The  other  two  bands  are  new.  We  find  that  the 

A 

waves  in  band  II  are  qualitatively  very  similar  to  the  Alfven  waves,  while 
the  waves  in  band  III  are  very  different. 

Finally,  we  estimate  the  total  radiated  power  for  a  "passive" 
condxting  sphere,  and  also  for  a  conducting  cylinder  oriented  such  that 
the  axis  of  the  cylinder  is  parallel  to  the  "motional"  electric  field.  In 
both  cases,  we  assume  that  the  source  current  density  is  uniform  everywhere 
inside  the  conducting  body,  and  that  it  is  in  the  direction  of  the 
"motional"  electric  field.  The  amplitude  of  the  source  current  can  then 
be  computed  using  Poynting's  Theorem  and  Ohm's  law.  Table  I  sumnarizes  the 
results  in  several  limits  where  simple  analytic  expressions  are  obtainable. 
In  particular,  for  the  case  of  a  wire  10km  long  and  less  than  Imm  in 
radius,  we  find  the  radiated  power  to  be  about  100  watts. 
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SECTION  1.  INTRODUCTION 


Phase  and  Group  Velocities  of  Waves  in  Cold  Plasma 


This  report  presents  results  of  a  comprehensive  analytical  and  numerical 
study  of  the  dispersion  relation  for  a  "cold"  plasma.  Although  the  subject 
matter  enjoys  a  time-honored  familiarity  for  anyone  interested  in  plasma 
waves,  it  appears  that  a  unified  presentation  of  phase  and  group  velocities 
for  all  frequencies  and  for  all  directions  of  wave  propagation  is  not  avail¬ 
able  in  the  literature.  Researchers  in  various  fields  dealing  with  plasmas 
(e.g.  the  earth's  magnetosphere,  the  interplanetary  medium,  the  solar  corona, 
etc.)  often  need  to  know  at  a  glance  the  magnitude  of  the  phase  or  the  group 
velocity  for  some  preferred  wave  frequency  and  d,*ection  relative  to  the  local 
magnetic  field.  The  most  convenient  way  to  supply  such  information  is  by 
means  of  polar  diagrams  (sometimes  also  referred  to  as  "Friedrichs'  diagrams") 
and  this  is  the  form  in  which  we  have  chosen  to  present  the  results  of  our 
calculations.  The  literature  is  replete  with  special  examples  of  such  dia¬ 
grams,  so  that  we  need  not  dwell  upon  the  pedagogical  details  of  their  con¬ 
struction.  Nonetheless,  for  the  convenience  of  the  reader,  we  do  summarize 
the  basic  equations  needed  for  the  derivation  of  the  dispersion  relation  (see 
Section  2).  In  Section  3  we  outline  the  mathematical  formalism  underlying  the 
concept  of  the  group  velocity.  Section  4  is  devoted  to  the  polar  diagrams  of 
both  the  phase  and  the  group  velocities.  A  brief  discussion  of  some  salient 
morphological  features  of  the  emerging  patterns  is  also  included.  Section  4 
contains,  mostly  in  graphical  form,  additional  information  that  elucidates  the 
causal  relations  between  the  phase  velocity  graphs  and  those  for  the  group 
velocities.  Section  5  addresses  in  .more  detail  the  dependence  of  our  results 
on  the  relative  magnitude  of  the  dimensionless  parameters  entering  into  calcu- 
lat'ons.  In  particular,  the  morphological  changes  of  the  polar  diagrams 
resulting  from  the  increase  of  the  magnetic  field  strength  are  presented  for 
various  cases  of  special  interest  in  astropnysics  and  solar  physics. 
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SECTION  2.  THE  DISPERSION  RELATION 


Phase  and  Group  Velocities  of  Waves  in  Cold  Plasma 


i 


Following  standard  procedures  (see,  e.g.,  Stix  1962)  we  shall^  assume  that 
the  plasma  is  in  the  prooer  frame  of  reference  (plasma  velocity  Vp  and  static 
electric  field  tg  a^-e  zero)  and  that  the  mean  velocity,  V|^,  of  the  particle  of 
type  "k.”,  the  electric  field  f,  the  current  density  j,  etc.  are  infinitesimal 
perturbations.  To  derive  the  dispersion  relation,  we  begin  with  the  linear¬ 
ized  equation  for  conservation  of  momentum  for  the  species  of  the  form; 


3V^/at  =  Z. a  (: 


\  X  Bg/C) 


(1) 


of  the  k' 


(1/-,,  w  «  V, 


where  Bg  is  the  external,  magnetic  field,  m|^  and  Z|^e  are  the  mass  and  charge 

particle,  respectively.  We  have  neglected  the  effect  of  pressure 
§o/c);  in  other  words  we  restrict  ourselves  to  the  cold 
plasma  aoproximation .  The  gravity  effect  is  unimportant  in  most  applications 
(notably  this  is  the  case  in  earth's  magnetosphere).  Solving  Eg.  (1)  for  V. 
allows  us  to  construct  the  expression  for  the  mobility  tensor  M,  which  we 
define  in  dimensionless  units  for  each  type  of  particle 


\  =  (c/8j) \  .  I 
where  M|^  is  given  by 


(2) 


1 

2  ^0)  +  £ 


k'-k 


k'‘k 


^k^-k 


^k^k 


“  *  "k^k 


2  u  +  CO  -  £^a^ 


•).  0 


-  I  .  ^k-k  ^  1  r 

2  CO  +  o  -  '  2 


^k“k 


^k‘\ 


r.  .  ^k^k 


).  0 


(3) 


gives  the  sign  of  the  charge,  and  *  I  (Zj^e  Bg)  /  (mj^c)!  is  the  cyclotron 
freouency  for  the  k^^  particle. 

Having  obtained  the  expression  for  the  mobility  tensor,  we  can  calculate 
the'  dielectric  tensor  K  defined  by 
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iT  •  E  =  E  +  (4ni/u))  j 


(4) 


where  the  current  dens’ity  is  given  by 
^  =  J  "k  ^k®  \ 

with  n|^  being  the  number  density.  In  terms  of  the  dimensionless  mobility  ten¬ 
sor,  Eq.  (2),  (3)  and  (4)  may  be  combined  to  yield 

T=T+  i  I  /  (u  Q,^)  (6) 

where  1  is  the  unity  matrix  and  the  plasma  frequency  of  the  k^^  species  is 
given  by 

"k  "  ^k^ 

To  obtain  the  dispersion  relation  we  calculate  the  electric  field  E  using 
Maxwell's  equations.  We  have 

a^E/3t^  r  7  *  (7  X  E)  =  -  4n  aj/at  .  (7) 


A'ter  Fourier  transformation  in  time  and  space  (assume  homogeneous  and  infi¬ 
nite  plasma),  Eq.  (7)  becomes 


k  * 


(k  X  E)  J/c^  T*  E  =  0 


(8) 


where  k  is  the  propagation  vector.  Equation  (8)  is  of  the  form  Ej  =  0, 
where  the  matrix  A.j  is  given  by 


2  2  2 
S-n  cos  e  ,  -ID  ,  n  cose  sine 


iD  ,  S-n  ,  0 


2  2  2 
n  cose  sine  ,0  ,  P-n  sin  e 


(9) 


The  symbols  S,  R,  L,  D,  P,  and  n  are  defined  in  a  conventional  way: 
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S  5  1/2  (R  ^  L);  0  =  1/2  (R  -  L) 

R  =  1  ~  Z(!Ij^^/w^)  u/(w  Ej^W)^) 

(10) 

L  =  1  '  E(ll,^^/u‘')  u/(u  - 

(11) 

p  —  ■<  «-’t2/2 

r-i'iil.  /u 
k  ** 

(12) 

-T  '  “  • 

■  . 

n  =  KC/id 

(13) 

two-..omoone''t  o’asma,  with  electrons  as  negative 
2)  reauce  to 

-ons,  Eqs.  (10),  (11) 

R  -  1  ~  /(w  C.j)(u  ~  Qg) 

(lA) 

L  =  1  -  Up^/(^  -  -.)(u.  *  Qg) 

(15) 

...  t 

P  =  1  - 

(16) 

the  plasma  frequency,  Up,  is  defined  as 

^  _  ,,  2  ^  „  2.  1/2 

Wp  -  (n,  *  ^ 

(17) 

« 

the  condition  for  a  nontrivial  solution  of  Eo.  (8)  we  obtain  the  familiar 
d^soersion  relation  for  the  phase  velocity  (for  references  see  Stix  (1962)). 
Wor<ing  out  the  algebra  of  det  (A.j)=0  we  find 


[(x-o)cos^9  -  x-o]  r  C«^c3[xsin^0  -  t/;cos^e]  =  0 


(18) 


where  U  =  w/k  is  the  phase  ve  ocity,  =  Bg/(4llp)^^^  is  the  Alfven  velocity 
for  the  case  in  which  the  magnetic  energy  is  substantially  less  than  the  rest 
energy,  and  9  is  the  angle  between  ^g  and  k.  The  dimensionless  quantities  x, 
a  and  ♦  are  defined  by 


X  =  -w^/  -  Wp^)) 

0  =  -  ((w^  -  ia,H^)  -  Wlc^) 


"UC 


(19) 

(20) 
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-  ^  c; 


-..I 


1 

J 


1 


(21) 


'P  =  ((J  -  -  “LH^)  ■  «Ur^)) 


wnere  u 


LH’  '^UH 


and  uyQ  are  the  lower  hyorid,  the  upper  hybrid  ana  the 


U.1  I 

lower  and  upper  cutoff  frequencies,  respectively.  They  are  aefined  as 


2 

“LH 


Ti-n  +  -'I- 
1  e  1 


2 

i_) 
2 
T  . 

1 


(22) 


2 

“LH 


(23) 


is  the  ratio  of  to  the  speed  of  light.  Note  that  x,  a  and  are  func¬ 
tions  of  ui  and  the  characteristic  plasma  frequencies,  but  do  not  depend  on  9. 
The  solution  of  Eq .  (18)  yields  for  U^/C^: 


=  i  3  ^[(x-o)  cos^e  -  x-o  t  Y(x'o)^  sin^e  +  4a(a  +  i)<)  cos^e]  (25) 
C?  2  A 

Equation  (25)  shows  that  there  are  only  two  modes  in  the  cold  plasma 
apo-ox’mation .  We  shall  call  them,  respectively,  the  "plus"  and  the  "minus" 
mode.  In  tne  literature  these  two  modes  go  by  a  variety  of  different  names 
qeoending  on  the  frequency  regime  one  happens  to  be  interested  in.  The  "plus" 
mode  is  variously  referred  to  as  the: 

*ast  mode  (0  <  ui  <  Q. ) 

lower  hybrid  mode  (il-  <  w  <  wy^) 

whistler  mode  (wy^  <  w  <  n^) 

electron  cyclotron  mode  (w  <  C^) 
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ordinary  mode  (u  > 

as  the  frequency  regime  is  changed;  similarly,  the  "minus"  mode  is  referred  to 
as  the; 

Alfven  mode  (0  <  u  < 

ion  cyclotron  mode  or  ion  whistler  mode  (w  <  fl-) 

upper  hybrid  mode  or  extraordi nary  mode  (w^  <  u  < 

In  order  to  unify  our  presentation,  and  simplify  the  nomenclature,  we  shall 
always  refer  to  these  modes  as  "plus"  and  "minus"  modes. 

The  dispersion  relation  (Eq.  18)  gives  zero  or  infinity  for  the  phase 
velocity  when  the  dimensionless  parameters  go  to  a  certain  value.  Stix  (1962) 
terms  the  case  U  0  a  resonance,  and  the  case  U  a  cutoff.  Resonances 
occur  when  the  following  criterion  is  satisfied: 


,2  ^  2..  2  „  2w  2  -  2s 

2  (u  -  )i^  -  ilg  )(u  ~  > 

tan  e  ■  - - j-j - J--2 - j- 

(i) 


(26) 


After  some  algebra  Eq .  (26)  becomes 


-  J  +  (1  -  COS^9)u.f(^]  + 


2.^  2 
LH- 


?  722  222-1 

+  u  [(1  -  COS^9)  (UyH  +  sip  Slg  COS  9]  - 


(27) 


where 


0  P  ’  n  ^ 
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The  square  roots  of  Eq.  (27)  are  easy  to  find  if  we  assume  that  w  <  or  w  < 
Solving  Eq.  (27)  for  these  conditions,  we  obtain  the  general  formulas 
for  lower  middle  (n^^j^)  and  upper  (!12y[^)  resonance; 


where 


^  2^  2  2  2, 

Q.  a  u)  cos  9 

1  e  p 

"X7T^~7T27T 


[(1  -  COS^9)u)y^  a)y^  +  COS^e] 


(28a) 


2  2 

^(l  *  cos  9)(Uy^ 


'MR 


2  2  2  2 
u  ■*■  COS  9]  + 

P  8 


2[a)J^  +  (1  -  COS^9)  ui^LH] 


(28b) 


A  =  [(1  -  COS^9)ouy^  a)y^  +  COS^9]^ 

-  4[u)yj^  +  (1  -  COS^9)uL^]a.^Qg^  u)p^  COS^9]  , 


1  2 
■  2  ‘^UH 


[1 


1  - 


cos^e  , 

4 -  J- 

“UH 


(28c) 


From  Equations  (28a),  (28b)  and  (28c),  we  can  obtain  the  resonant  frequencies 
for  tne  part;cu1ar  physical  situations. 


The  cutoffs  occur  when  the  following  conditions  are  fulfilled: 


or 


or 


4,  -  - 


Table  1  shows  the  dependence  of  tne  resonances  and  cutoffs  on  the  angle  9, 
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The  ordering  pre'.ented  in  this  sketch  (see  Table  2)  is  typical  for  weak 

2  2  2 

magnetic  fields  (w  p  »  2^  or  <  0.001),  and  change  significantly  when  the 
magnetic  field  increases  (see  Section  6  for  a  more  detailed  discussion). 
Table  2  presents  values  of  the  resonances  and  cutoffs  (8  =  0)  for  earth's  mag¬ 
netosphere  =  0.001  and  electron  density  n^  =  3.6x10^). 

TABLE  2 

The  resonances  and  cutoffs  for  the  earth's  magnetosphere 

5  -3 

with  8^  =  0.001,  B  =  0.336  and  n^  =  3.6x10  cm  . 


“lh 

^e 

“LC 

“p 

“UH 

“uc 

1.96x10^ 

3.33x10^ 

5.80x10® 

3.11x10^ 

3.39x10^ 

3.44x10^ 

3.69x10^ 

Phase  and  Group  Velocities  of  Waves  in  Cold  Plasma 


AO 


SECTION  3.  GROUP  VELOCITY 


The  group  velocity  of  a  wave,  U^,  is  defined  in  general  by  equation 


Ug  = 


(29) 


In  dispersive  and  anisotropic  media  the  direction  and  the  magnitude  of  Ug  are 
different  from  those  of  U.  The  sketch  below  defines  the  pertinent  symbols 
relating  Ug  to  U. 


^  A,  ^ 

Since  w  may  be  locked  upon  as  a  function  of  k.j  and  (k  =  k||k||  *  k_i_kj_)  Eq. 
(29)  imolies  that 


lr„  ■''11  • 


(3o: 


For  the  explicit  evaluation  of  3w/3kj_  and  3«/3k|j  we  need  the  dispersion  rela¬ 
tion 


F(a)*  k||  y  k^)  -  0 


(31 : 


obtainable  from  Eq.  (18).  We  then  have 
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7*  'W  1^  '^'B* 


+ 


If. 

3kf| 


dk„ 


(32) 


and 


U  =  -  ( 

n  ' 


+  If.  k  )  _L 


"aFii  II  akj. 


(33) 


where  F  =  3F/3w.  For  the  modes  considered  here  we  find  from  Eqs.  (32)  and 
w 

(33)  the  following  expression  for  the  group  velocity: 


r2 

Cac 


[cj  in  cos^e  -  X  sin^9)-U^]  + 


k  f 

*  [C,^*  -  f  -  J  [Ca^o  *  T  i 


(34) 


where 


F  H  k"^  F  =  2U^  +  UC.^[2a  cos^e  +  (x  -i-  j)sin^9  + 

U  uj  A 


c 

*  j  {o  *  <3q) {]  *  cos^e  +  (x  sin^e  -ih  cos^e))  + 


+  (1  -  xSa  )(1  +  a  )  X  sin  9 


2 

-  - J  ('1'^  ■*■  ('g)<7  COS  0] 


(35) 
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^  (Zo)  +  fig  -  n,-) 

^  [n.  Hg  -  -  n.)(y  +  Qg)]^ 


'^B  ■  2 


^  wfl.ng  (2a.  -  Jig  + 


[n,-  flg  -  8^^(a.  -  n^)(aj  -  Sig)]^ 


01^  n.n 
1  e 


(a."  -  a.")  (a."  -  q/) 


(ai-Jl.I5g)  (2a)^  -  fig^  -  ) 

(J  -  n.^){i/  -  ag^) 


Introducing  the  polar  angle  9g  shown  in  the  sketch  above,  we  can  calculate, 
according  to  Eq.  (34),  the  components  of  Ug  parallel  and  perpendicular  to  Bq. 
Then  we  have 


C  ^0 

U„  z  U  COS0  =  -p —  [C.^(ij;(l  +■  cos^e)  -  X  sin^e)-2U^]  cose 

3  3  '  11  ^ 


U,  =  sine^  =  4-^  CC.^(-4/  cos^e  -  x(l  +  sin^9))-U^(l  -  J  )]  sit 

i  g  g  M  o 


The  above  equations  yield  then,  with  the  help  of  Eq.  (18), 


-  U^C.^  (x  -  <^)  cos^e  -  C.^  ox 

tan  9  =  -T - 5 — 5 - ; - 2 - 5 —  ® 

9  (x  -  o)  sin^e  +  ax” 
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and 

Ug  =  -V  [U®  +  j  (X  -  a)[^  u^c/  (X  -  a)  +  (,<  +  '^)]  sin^20  + 

u 


+  C.^a  [ZU^  (x  +  !(;)  +  C.*^c  (ij;^  -  X^)  C0S^9  + 


(38) 


+•  C^^<^X  -  2U^]]. 


Since  the  general  equations  for  the  group  velocity  are  cumbersome,  it  is 
instructive  to  derive  simplified  expressions  for  and  in  the  MHO  limit 
(u  <<  S2^).  Thus  we  get  x  "*■  0 ,  o  ~  -I,  ♦  ~  1,  (o^  +  Og)  0,  +  i/»g)  -*  0  and 


which  yields 


9 


U 


g 


2 


tan  9g 


tan  9 


for  the  fast  mode 


(39) 


and 


Ug^  =  ®g  ~  ^  Alfven  mode. 
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The  solutions  of  the  dispersion  relation  for  the  phase  and  group  velocity  • 
(Eqs.  25,  37  and  38)  are  plotted  as  the  po'ar  diagrams.  For  numerical  calcu¬ 
lations  it  is  convenient  to  redefine  the  dimensionless  parameteri  that  were 

introduced  in  Eqs.  (25),  (37)  and  (38).  After  defining  the  mass  ratio  as  ?  = 
1/2 

(m^/m^)  ,  we  have 


X  E 


(40) 


1 

a  =7 


1 


(f  -  i;)(l  +  fc)  ' 


1 


c 

(f  +  Jcr-TF) 


(41) 


Ip  =  -s — - 

;2(f2  .  1) 


(f^  -  C^)(l  - 


(42) 


and 


=  1 


1  -  -S- 


A=2  (f  -  ;)(1  >  f^T  [7  -  sf  (f  -  c)(l  +  fc)] 
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(46) 


=  (1  - 


nl) 


2f^ 


{r  -  1) 


-  2r 


where  f  =  Now  Eqs.  (25).  (37),  (38)  and  (40)*(46)  represent  a 

I  C 

convenient  form  for  calculations  because  U,  and  0^  are  functions  of  u  and 
two  physical  parameters  and  Thus,  holding  B^^  and  C  fixed  we  construct  a 
family  of  polar  diagrams  for  U  and  U  for  various  frequencies.  We  plot  the 
polar  diagrams  using  different  normalization  of  f  for  different  regimes  of 
frequency: 

I  and  II:  f  =  f-  C  =  («-/«,•  )4 

III;  fa  s  (w/wlh)  wr  =  («/««i_h)  ((^i^  *  H.^)  /  (Q^^  +  + 

IV:  f  5  fg  s  («/fig)c‘^ 

V  and  VI;  f  5  fpfl^"^  5  (y/u^) 

The  results  presented  in  this  section  are  typical  for  the  earth's  magneto- 
spnere,  where  c  =  1/171.4  and  6^  =  lO”^.  Figures  4.1  through  4.5  show  the 
polar  diagrams  for  the  phase  velocity  of  the  "plus"  mode  in  “.he  different  fre- 
quency  regions.  The  diagrams  for  the  group  velocity  of  th»  "plus"  mode  are 
presented  on  Figures  4.6  through  4.10.  The  next  four  figures  show  the  dia¬ 
grams  for  U  and  of  the  "minus"  mode.  As  one  sees  from  Figures  4.1  through 
4.14,  there  are  several  general  trends  in  the  behavious  of  the  modes.  The 
polar  diagrams  for  the  "plus"  mode  present  the  "evolution"  of  loci  of  the 
phase  velocity,  from  a  sphere  (Fig.  4.1)  through  two  distinct  spheres  that  are 
surfaces  of  revolution  about  the  symmetry  axis  (Fig.  4.2);  these  evolve  even¬ 
tually  to  small  aspect  ratio  lobes  (Figs.  4.3a  and  4.3b)  which  are  narrowly 
confined  in.  The  symmetry  axis  is  defined  by  direction  of  the  magnetic  field 
that  is  along  the  abscissa  axis.  For  frequencies  just  above  W|_0,  a  sphere 
appears  again  (Fig.  4.4)  and  becomes  slightly  deformed  in  the  regime  VI  (Fig. 
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4.5).  The  loci  of  the  group  velocity  evoUe  from  a  sphere  (Fig.  4.6)  through 
a  "pre-whistler"  shape  (Fig.  4.7)  to  the  rharacteri stic  whistler  cusp  (Fig. 
4.8a).  However,  the  behavior  of  the  loci  'n  the  regime  is  more  complex  than 
in  the  other  frequency  regions.  The  group  velocity  locus  can  be  plotted  as 
one  topologically  complex  curve,  but  for  prjsentation  it  is  convenient  to  dis¬ 
tinguish  two  different  shapes;  the  whistlir  cusp  (its  evolution  is  shown  on 
Fig.  4.8a)  and  the  "expanding  lobes"  (Fig.  4.8b).  When  the  frequency 
increases  the  whistler  cusp  becomes  deformed  and  disappears  completely  at  w  = 
i2g/2;  as  Fig.  4.8c  shows  for  higher  frequenzies  the  lobe  dominates  (region  III 
is  discussed  later  in  this  section).  The  g-oup  velocity  loci  in  regimes  V  and 
VI  evolve  from  a  sphe'-e  (Fig.  4.9)  to  a  tiny  shape  limited  to  the  direction 
perpendicular  to  ?  (Fig.  4.10). 


The  next  four  figures  demonstrate  the  behavior  of  the  "minus"  mode.  F’g- 
ures  4.11  and  4.12  show  that  the  phase  velocity  does  not  change  its  character¬ 
istic  shape,  but,  that  during  the  evolution  the  loci  shaoe  becomes  smaller  and 
disappears  at  the  critical  frequenc’es  and  «y|^.  The  group  velocity  loci 
(Figures  4.13  and  4.14)  show  the  ion  whistler  cusp  in  the  regime  I  and  the 
cone  shape  (regime  VI)  that  is  deformed  when  the  frequency  approaches  wy|^.  The 
angle  9^  reaches  the  maximum  value  at  12°  (for  f^  -  0.85)  in  regime  I  and  the 
minimum  value  at  12°  (for  f^  1.005)  in  regime  VI. 


The  "plus"  mode  is  the  only  one  that  can  exist  in  the  region  III;  'n  this 
region  the  group  velocity  loci  have  topologically  complex  structure  that 
strongly  depends  on  frequency.  The  series  of  Figures  4.15  through  4.25  pres¬ 
ent  the  changes  of  U  and  (plotted  togetier)  for  increasing  frequency.  One 
sees  the  characteristic  whistler  cusp  :  1);  where  the  line  that  is  plotted 
along  B  (and  ties  the  origin  to  the  ends  ot  the  cusps)  marks  the  iniiial  stage 
of  formation  of  the  lobes  (Figures  4.15  and  4.16).  Generally  this  "line"  is 
split  and  the  split  angle  is  as  small  als  0.017°  for  =  1.01,  When  the  fre¬ 
quency  increases,  the  split  angle  increases  also,  and  the  lobe  becomes  the 
significant  feature  on  the  diagrams  (F^g.  4.17  through  4.22)  instead  of  the 
whistler  cusp.  The  critical  point  occt'rs  for  u  =  S2./2  (or  f ,  m  =  86)  when  the 
whistler  cusp  disappears  (Fig.  4.23)  'i\no  expanding  lobe  can  be  seen  only 
(Figs.  4.24  and  4.25).  Therefore,  it  •.(;  convenient  to  divide  the  region  III 
into  two  subregions:  III4  (u  <  12^/2)' and  IIIB  (w  >  52^/2).  The  angle  9 
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reaches  the  maximum  at  76.5® 
however  the  maximum  angle  for 
gion  IIIA). 


for  fjj^  =  170  or 
the  whist’er  cusp 


fg  =  0.98  (subregion  IIIB) 
never  exceeds  19.29°  (subre 
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SECTION  5.  DEPENDENCE  OF  THE  DISPERSION  RELATIONS  ON  THE 


MAGNETIC  FIELD  STRENGTH  AND  DENSITY  OF  PLASMA 


and  Group  Velocities  of  Waves  1n  Cold  Plasma 


The  general  equations  for  the  phase  and  group  velocity  (Eqs.  25,  37  and  38) 
contain  only  two  dimensionless  parameters:  c  and  These  parameters  influ¬ 
ence  the  results  by  changing  the  values  of  the  characteristic  frequencies  and* 
thus  changing  the  widths  of  the  regimes.  Generally,  C  affects  all  the  charac¬ 
teristic  frequencies  except  Sg  and  Wp  and  consequently  the  change  of  the  fre¬ 
quency  regimes.  However,  the  change  of  not  only  implies  the  change  of  the 
width  of  the  regimes  but  can  also  lead  to  a  redefinition  of  the  upper  and 
lower  limits  of  these  regimes.  The  cutoff  frequencies  anc  are  the 
most  important  parameters  for  changing  the  width  of  the  regimes.  Table  3 
shows  the  dependence  of  on  0^  and  c  and  the  "vanishing"  regime  in  the  var¬ 
ious  cases. 

TABLE  3 


Dependence  of  the  cutoff  on  the  changes  of  8^^ 


4 

“LC 

2 

“P 

Vanishing 

regime 

-■“P 

. 

» 

V 

^MR  “  ^e 

ZQ  * 
e 

IV 

2 

3 

'“LH 

iJg  V  .j  f2g 

f 

III 

1.0 

\r  * 

II 

2.0 

0 

7  ^*i“e 

I 
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In  Section  2  we  defined  the  nine  frequency  regimes  for  the  particular  phys¬ 
ical  conditions  typical  for  the  earth's  magnetosphere.  Now,  we  are  interested 
in  how  the  ordering  and  the  classification  will  be  altered  when  8^  increases. 
Table  4  presents  a  sequence  of  inequalities  for  and  shows  the  changes  of 
the  bounds  of  the  regimes.  As  one  sees  from  Table  4  a  significant  reordering 
of  the  regimes  is  necessary  as  8  increases. 

TABLE  4 

Dependence  on  the  bounds  of  the  frequency  regime  on  8^. 


1 

Regime 

3^  <  0.001 

0.001  <  <  0.01 

0.01  <  sj  <  0.1 

2 

0.1  <  <  1.0 

I 

0,  n^R 

0, 

0,  nj_R 

0, 

II 

\r*  “lh 

\r*“lh 

\r*  "lc 

Vanishes 

III 

“LH’  ^ 

"LH*“LC 

Vanishes 

Vanishes 

i 

i 

I 

SlR‘  “LC 

Vanishes 

Vanishes 

Vanishes 

i 

! 

i  V 

1 

“LC’  "p 

“LC’  ‘"p 

“LC’  “p 

“LC’  \r 

VI 

“p’ 

“p’  ^MR 

“p»  ^r 

“p*  ^ 

Let  us  consicer  the  results  of  calculations  that  were  made  for  varying 
between  0.0005  and  0.05  and  for  t  *  1/171.4  (for  oxygen).  The  polar  diagrams 
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for  the  group  velocity  are  presented  on  Figs.  5.1  and  5.2,  and  show  the 
independence  of  the  shape  of  the  loci  on  as  long  as  the  regime  is  not  dis¬ 
turbed  by  Figure  5.2  presents  the  disturbance  of  regime  III  by  that 
becomes  the  new  "bound"  for  this  regime.  For  frequencies  higher  than  but 
less  than  Qg  the  behavior  of  the  group  velocity  is  typical  for  regime  V  (see 
Section  4). 

Now,  we  consider  the  results  for  two  different  mass  ratios:  =  1/171.4 

(for  oxygen)  and  =  1/42.85  (for  hydrogen)  and  for  0^  =  0.1,  0.03,  0.01, 

0.005  and  0.001  whe-e  the  changes  of  and  B  are  treated  independently.  The 

values  of  6^  that  were  chosen  can  be  identified  with  different  astrophysical 

objects  such  as  the  atmospheres  of  magnetic  white  dwarfs,  the  atmospheres  of 

magnetic  (Ap)  stars,  and  the  solar  atmosphere  (e.g.,  chromosphere  and  coronae 

above  sunspots).  The  results  are  presented  on  0_  versus  9  diagrams  that  give 

direct  information  about  the  connections  between  the  direction  of  k  and  U^. 

0 

The  series  of  the  ciagrams  (Figs.  5.3  through  5.6  for  hydrogen  and  Figs.  5.7 
through  5.10  for  o<ygen)  are  plotted  for  the  two  constant  frequencies  u  = 
5x10^  and  w  =  5x10^.  In  addition,  the  calculations  were  carried  out  for  the 
given  0^  but  with  fonstant  n^  (=  1 .Oxl.O^^cm”^)  or  constant  B  (=  5  G).  The 
mode  frequency  was  leld  fixed  and  the  values  of  B  and  n^  were  chosen  for  con¬ 
venience  of  presentation  for  the  various  frequency  regimes.  The  9g  versus  8 
diagrams  (Figs.  5.3  through  5.6  and  5.7  through  5.10)  present  changes  of  the 
behavior  of  the  group  velocity  for  different  0^;  a  comoarison  of  Figs.  5.3 
with  5.7  and  5.4  with  5.8  shows  the  results  for  the  different  c.  Several  fea¬ 
tures  of  these  diagrams  merit  special  emphasis: 

1.  the  "formation"  of  the  whistler  mode  from  the  almost  ideal  MUD  fast  mode 
for  decreasing  0^^  (Figs.  5.4,  5.7  and  5.8),  or  the  reverse  process  for 
increasing  0^; 

2.  a  similar  evolution  can  be  seen  for  decreasing  c  (compare  Figs.  5.3  with 
5.7  or  5.4  with  5.8) ; 

3.  the  absence  of  one  or  two  curves  in  Figs.  5.5,  5.6,  and  5.9  is  caused  by 
the  independence  of  the  shape  of  the  loci  on  the  specific  values  of  0^; 
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4.  the  crossing  of  the  bound  between  regimes  III  and  V  is  shown  on  Flg.  5.6 
and  5.9;  Fig.  5.10  presents  the  crossing  of  the  bounds  between  regimes 
III-V  and  V-VI  (regime' IV  does  not  exist  for  >  0.001). 
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In  this  report  we  derived  the  general  dispersion  relations  (Eqs.  25,  37  and 
38)  for  the  phase  and  group  velocity  in  the  “cold"  plasma  approximation.  The 
dispersion  relations  were  obtained  as  a  function  of  the  physical  parameters: 
C,  5^,  frequency  and  the  direction  of  wave  propagation.  In  this  approxima¬ 
tion,  there  are  only  two  clearly  d’stinct  modes:  "plus"  and  "minus".  Using 
the  general  definitions  of  the  resonant  (Eqs.  28a,  b,  c)  and  cutoff  (Eqs.  22, 
23  and  Z4)  frequencies,  we  classified  the  various  frequency  regimes,  and 
showed  in  which  regimes  the  modes  become  evanescent.  The  general  definitions 
of  the  resonant  and  cutoff  frequencies  allowed  us  to  show  their  dependence  on 
the  direction  of  propagation  (Table  1)  and  on  the  physical  parameters  (Table 

4). 


The  polar  and  0^  versus  9  diagrams  that  are  presented  in  the  report  show 
the  dist’nct  behavior  of  the  "plus"  and  "minus"  modes  in  different  frequency 
regimes.  Each  particular  frequency  domain  is  characterized  by  a  typical  shape 
of  the  phase  or  group  velocity  loci.  This  shape  is  independent  of  the  physi¬ 
cal  parameters  c  and  (Eig.  26)  that  change  the  width  of  the  regimes  in  fre¬ 
quency  space,  the  identity  of  the  regimes,  as  well  as  of  the  value  of  the 
phase  or  group  velocity. 

A  particular  physical  condition  (B^^  *  0.001)  gives  the  maximum  number  (nine 
according  to  our  definition)  of  distinct  regimes  of  behavior.  The  "plus"  mode 
can  propagate  in  eight  of  these  frequency  domains,  but  the  propagation  of  the 
"minus"  mode  is  admissible  in  only  three  of  them.  There  is  only  one  regime 
(IV)  where  both  modes  are  evanescent.  Either  an  increase  or  decrease  in  B^ 
causes  a  decrease  in  the  number  of  distinct  regimes  of  behavior. 

In  region  II,  we  showed  (Fig.  4.7)  that  the  "plus"  mode  exhibits  a  degener¬ 
acy  in  group  velocity  (i.e.,  there  exist  three  distinct  group  velocities  for 
one  very  specific  angle  9^);  this  degeneracy  never  occurs  in  the  direction  of 
the  magnetic  field.  One  can  also  see  the  degeneracy  for  one  angle  9^  in  Figs. 
4.20  and  4.21.  However,  in  these  latter  cases,  the  physical  picture  is  more 
complicated,  and  will  be  discussed  in  a  separate  paper.  Generally,  the  "plus" 
mode  never  exhibits  the  degeneracy  for  one  0^  except  in  some  portion  of 
Region  II. 


Phase  and  Group  Velocities  of  Waves  in  Cold  Plasma 


Several  characteristic  angles  0^  were  found  for  the  modes.  The  transport 
of  energy  by  the  "plus"  mode  is  restrictd  to  angles  19.29®  and  76.5®  for  sub- 
regimes  III  A  and  III  B,  respectively.  The  minus  mode  can  transport  energy  at 
less  than  the  angle  0^  =  12®  in  Regime  I  and  only  above  the  angle  0^  =  12®  in 
Regime  VI.  The  coincidence  will  be  also  discussed  in  a  separate  paper. 
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FICITRE  CAPTIONS 


Figure  4.1 


The  "plus"  mode  in  the  earth's  magnetosphere.  We  show  the  polar 
diagrams  for  the  phase  velocity  in  Regime  I  and  for  frequencies 
f^.  =  0.05,  0.35,  "O.S,  0.75,  0.85  and  1.0  (correspond  with  the 
curves  from  the  internal  to  the  external  one,  respectively).  The 
magnetic  field  is  along  the  abscissa  axis. 


Figure  4.2  The  same  as  Fig.  4.1  but  for  Regime  II  and  for  frequencies  f^-  = 
1.1,  3.5,  8.5,  17.0,  25.0,  33.0  and  41.0. 

Figure  4.3a  The  same  as  Fig.  4.1  but  for  subregime  III  A  and  for  frequencies 
f^H  =  1-1,  2.0,  5.0,  10.0,  25.0.  50.0  and  86.0. 

Figure  4.3b  The  same  as  Fig.  4.1  but  for  subregime  III  B  and  for  frequencies 
f^^  =  170.0,  150.0,  120.0,  100.0  and  86.0. 

Figure  4.4  The  same  as  Fig.  4.1  but  for  Regime  V  and  for  frequencies  f  = 
0.96,  0.95,  0.94,  0.93  and  0.92. 

Figure  4.5  The  same  as  Fig.  4.1  but  for  regime  VI  and  for  frequencies  f^j|^ 
=  0.995,  0.994,  0.993.  0.992,  0.991  and  0.990. 

Figure  4.6  The  plus  mode  is  the  earth's  magnetosphere.  We  show  the  polar 
diagrams  for  the  group  velocity  in  Regime  I  and  for  frequencies 
f.j  =  0.05,  0.35,  0.5,  0.75,  0.85  and  1.0  (correspond  with  the 
curves  from  the  internal  to  the  external  one,  respectively). 
The  magnetic  field  is  along  the  abscissa  axis. 

Figure  4.7  The  same  as  F-g.  4.6  but  for  Regime  II  and  for  frequencies  f.j  = 
1.10,  3.5,  8.5,  17.0,  25.0,  33.0  and  41.0. 

Figure  4.8a  The  same  as  F  g.  4.6  but  for  subregime  III  A  and  for  frequencies 
^LH  ”  l^-0.‘'25.0  and  50.0.  The  evolution  of  the  whis¬ 

tler  cusp  only  is  presented. 

Figure  4.8b  The  same  as  F’g,  4.6  but-  for  subregime  III  A  and  for  frequencies 
flH  “  50.0.  The  splitting  of  the 

"line"  for  increasing  frequency  is  shown. 
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Figure  4.8c  The  same  as  Fig.  4.6  but  for  subregime  III  B  and  for  frequencies 
f^H  =  170.0,  150.0,  120.0,  100.0  and  86.0. 

Figure  4.9  The  same  as  Fig.  4.5  but  for  regime  V  and  for  frequencies  fy^  = 
0.96,  0.95,  0.94,  0.93  and  0.92. 

Figure.  4.10  The  same  as  Fig.  4.6  but  for  regime  VI  and  for  frequencies  fy^  = 
0.995,  0.994,  0.993,  0.992,  0.991  and  0.990. 

Figure  4.11  The  "minus"  mode  in  the  earth's  magnetosphere.  We  show  the 

polar  diagrams  for  the  phase  velocity  in  regime  I  and  for  fre¬ 
quencies  f^  =  0.99,  0.95,  0.85,  0.75,  0.50,  0.35  and  0.05 
(decreasing  from  the  internal  to  the  external  curves,  respec¬ 

tively).  The  magnetic  field  is  along  the  abscissa  axis. 

Figure  4.12  The  same  as  Fig.  4.11  but  for  regime  VI  and  for  frequencies  fy^ 
=  0.995,  0.994,  0.992,  0.990  and  0.899. 

Figure  4.13  The  "minus"  mode  in  the  earth's  magnetosphere.  We  show  the 

polar  diagrams  for  the  group  velocity  in  the  regime  I  and  for 
frequencies  f^  =  0.99,  0.95,  0.85,  0.75,  0.50,  0.35  and  0.05 
(decreasing  from  the  internal  to  the  external  curves,  respec¬ 

tively).  The  magnetic  field  is  along  the  abscissa  axis. 

Figure  4.14  The  same  as  Fig.  4.13  but  for  regime  VI  and  for  frequencies  fy^ 
=  0.995,  0.994,  0.992,  0.990  and  0.899. 

Figure  4.15  The  "plus"  mode  in  the  earth's  magnetosphere.  We  show  the  polar 
diagrams  for  the  phase  and  group  velocities  in  regime  III  and 
for  frequency  fy^  =  1.01.  The  magnetic  field  is  along  the 

abscissa  axis. 

Figure  4.16  The  same  as  Fig.  4.15  but  for  fyH=  1.1. 

Figure  4.17  The  same  as  Fig.  4.15  but  for  =  2.0. 
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Figure  4.18  The  same  as  Fig.  4.15  but  for  =  10.0. 

Figure  4.19  The  same  as  Fig.  4.15  but  for  =  15.0. 

Figure  4.20  The  same  as  Fig.  4.15  but  for  =  20.0. 

Figure  4.21  The  same  as  Fig.  4.15  but  for  f^^^  =  50.0. 

Figure  4.22  The  same  as  Fig.  4.15  but  for  fj^^  =  70.0. 

Figure  4.23  The  same  as  Fig.  4.15  but  for  fj^^  =  86.0. 

Figure  4.24  The  same  as  Fig.  4.15  but  for  f^^j^  =  120.0. 

Figure  4.25  The  same  as  Fig.  4.15  but  for  f^^j^  =  170.0. 

Figure  5.1  Propagation  of  the  "plus'*  mode  with  constant  frequency  u  = 
1330.0  s"^  (Regime  II)  through  the  earth's  magnetosphere  with  6^ 
=  0.005,  0.01,  0.015,  0.025,  0.035,  0.040  anr^  0.050  but  for  con¬ 
stant  the  magnetic  field,  8  =  0.33G.  Increasing  runs  from 
the  internal  to  the  external  curves,  respectively.  The  magnetic 
field  is  along  the  abscissa  axis. 

Figure  5.2  The  same  as  Fig.  5.1  but  for  w  =  8.5x10^  s  (Regime  III).  The 
change  of  the  shape  of  the  loci  during  the  propagation  is  caused 
by  the  crossing  of  the  bound  between  regimes  III  and  V  < 

S2e)- 

Figure  5.3  The  9^  versus  9  diagram  for  the  "plus"  mode  with  constant  fre¬ 
quency  w  =  5x10^  s~^.  The  curves  (from  the  bottom  of  figure) 

present  the  dependence  of  9„  on  9  for  the  values  of  B.  =  0.001, 

9  ^ 

0.005,  0.01,  0.03,  and  0.1  in  a  hydrogen  atmosphere.  The  number 
density  is  constant  and  equal  to  n^  =  1.0x10^*^  cm 

Figure  5.4  The  same  as  Fig.  5.3  but  for  w  =  5x10®  s 
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Figure  5.5  The  same  as  Fig.  5.3  but  for  B  =  5  G. 


Figure  5.6  The  same  as  Fig.  5.3  but  for  w  =  5x10^  s  ^  and  for  B  =  5  G. 

Figure  5.7  The  same  as  Fig.  5.3  but  for  an  oxygen  atmosphere. 

Figure  5.8  The  same  as  Fig.  5.3  but  for  u  =  5x10^  s‘^  and  for  an  oxygen 

atmosphere. 

Figure  5.9  The  same  as  Fig.  5.3  but  for  B  =  5  G  and  for  an  oxygen  atmos¬ 
phere. 

Figure  5.10  The  same  as  Fig.  5.3  but  for  u  =  5x10^  and  for  B  =  5  G  and  for 
oxygen  atmosphere. 
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